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AUTOMATED OPTIMIZATION OF MULTI- 
ELECTRODE PACING FOR CARDIAC 
RESYNCHRONIZATION 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] Pursuant to 35 U.S.C. § 1 19 (e) 5 this application claims priority to the 
filing dates of: United States Provisional Patent Application Serial No. 60/667,575 filed 
March 31, 2005; United States Provisional Patent Application Serial No. 60/667,529 
filed March 31, 2005; United States Provisional Patent Application Serial No. 
60/684,751 filed May 25, 2005; and United States Provisional Patent Application Serial 
No. 60/695,577 filed June 29, 2005; the disclosures of which are herein incorporated by 
reference. 

BACKGROUND 

Field of the Invention 

[0002] The present invention relates to monitoring and administering electric 
signals in living tissues. In particular, the present invention relates to a method and 
system for automatically optimizing multi-electrode pacing for cardiac resynchronization 
therapies. 

Related Art 

[0003] Cardiac resynchronization therapy (CRT) is an important new medical 
intervention for patients suffering from congestive heart failure. In congestive heart 
failure, symptoms develop due to the inability of the heart to function sufficiently well as 
a mechanical pump to supply blood to meet the body's physiological needs. Congestive 
heart failure is typically characterized by a gradual decline in cardiac function, 
punctuated by severe exacerbations leading eventually to death. Over five million 
patients in the United States are estimated to suffer from this malady. 
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[0004] The aim of cardiac resynchronization pacing is to induce the 
interventricular septum and the left ventricular free wall to contract at approximately the 
same time. Resynchronization therapy seeks to provide a contraction time sequence that 
can most effectively produce maximal cardiac output with minimal total energy 
5 expenditure by the heart. Generally, the optimal timing may be calculated by reference to 
hemodynamic parameters such as dP/dt, the first-order time derivative of the pressure 
waveform in the left ventricle. The dP/dt parameter is a well-documented proxy for left 
ventricular contractility. 

[0005] In current practice, external ultrasound measurements are used to 

1 0 calculate the dP/dt parameter. An external ultrasound is used to observe wall motion 
directly. A clinician uses the ultrasound system in a tissue Doppler mode, which 
provides a feature known as tissue Doppler imaging (TDI), to evaluate the time course of 
displacement of the septum relative to the left ventricle free wall. Ultrasonographic 
evaluation using TDI or a similar approach can be an important tool for qualifying 

15 patients for CRT. 

[0006] CRT as currently delivered is effective in about half to two-thirds of 
patients implanted with a resynchronization device. In approximately one third of these 
patients, CRT provides a two-class improvement in patient symptoms as measured based 
on the New York Heart Association (NYHA) fair level scale. In another one third, a one- 

20 class improvement in cardiovascular symptoms can be accomplished. In the remaining 
third, there is no improvement or, in a small minority, even a deterioration in cardiac 
performance. This third group of patients is referred to as non-responders. Furthermore, 
the NYHA one-class responders may experience only marginal or partial improvements, 
given the unfavorable results seen in the minority. 

25 [0007] The synchronization therapy, in order to be optimal, targets the point of 

maximal delay on a cardiac wall segment, and advances the timing thereof to synchronize 
contraction with an earlier contracting region of the heart, such as the septum. The 
current placement technique for CRT devices is mostly empirical. A physician typically 
cannulates a vein that appears to be in the region described by the literature as most 

30 effective. He can then position the pacing device and perform stimulation. The lack of 
extra cardiac stimulation, such as diaphragmatic pacing, is confirmed. The current 
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available techniques rarely provide sufficient time or means for optimizing the cardiac 
resynchronization process. 

[0008] In general, CRT optimization must be performed through a laborious 
manual method performed by an ultrasonographer evaluating cardiac wall motion at 
different lead positions and different interventricular delay (IVD) settings. IVD refers to 
the difference between the timing of pacing pulses sent to different locations, such as the 
right ventricle and left ventricle. Pace makers can vary the atrio-ventricular delay, which 
is the delay between stimulations sent to the atria and the ventricle or ventricles. These 
timing settings are important in addition to the location of the pacing electrode in the left 
ventricle. 

[0009] Some research efforts to assess cardiac motion through internal sensors 
have been made in the past. U.S. Patent No. 6,625,493, entitled "Orientation of Patient's 
Position Sensor Using External Field," to Kroll et. al, filed on Aug. 24, 2001, discloses a 
positional accelerometer for indicating a patient's vertical acceleration. U.S. Patent No. 
6,002,963, entitled "Multi-Axial Accelerometer-Based Sensor for an Implantable 
Medical Device and Method of Measuring Motion Measurements There from," to 
Mouchawar et. al, filed on Feb. 17, 1995, discloses cardiac wall motion detection using 
an accelerometer. U.S. Patent No. 5,991,661, entitled "System and Method for 
Measuring Cardiac Activity," to Park et. al, filed on Oct. 17, 1997, discloses the use of an 
accelerometer for rate-adaptive pacing. U.S. Patent No. 6,044,299, entitled "Implantable 
Medical Device Having an Accelerometer," to Nilsson, filed on Sept. 19, 1997, discloses 
an in-can accelerometer for providing rate control. 

[0010] In addition to incorporating sensors in a pacing lead, a pacing lead can 
accommodate multiple pacing electrodes to provide flexibility in selecting pacing 
locations. U.S. Patent 6,473,653, entitled "Selective Activation of Electrodes within an 
Implantable Lead," to Schallhorn et al., filed on Mar. 2, 2000, discloses an implantable 
multi-electrode lead adapted to allow selective activation of the included electrodes to 
electrically excite the tissue in the vicinities of the activated electrodes. The Schallhorn 
system requires a controller for each electrode, and separate conductors for sending 
analog excitation pulses configuration commands to these controllers. Consequently, 
these controller circuits and conductors make the Schallhorn system bulky and unsuitable 
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for many size-sensitive applications. U.S. Patent 5,593,430, entitled "Bus System for 
Interconnecting an Implantable Medical Device with a Plurality of Sensors," to Renger, 
filed on Jan. 27, 1995, discloses a two-conductor bus system for connecting physiologic 
sensors to a pacemaker. The two-conductor bus provides power to the sensors, and the 
sensors 5 output signals are modulated on the two wires. 

[0011] With the flexibility of a multi-electrode lead comes the difficulty in 
selecting the optimal electrode and timing configurations, because the number of possible 
configurations increases geometrically with the number of available electrodes. 
Currently, TDI is the common technique for characterizing dyssynchrony with accuracy. 
TDI measures myocardial movement throughout the cardiac cycle and can quantify the 
timing differences between peak tissue movements in multiple areas of the left ventricle, 
which suggest the presence of dyssynchrony. TDI can also be used to determine cardiac 
wall positions via external ultrasonography, typically for purposes of measuring valve 
function, cardiac output, or rarely, synchrony index. However, TDI is a complex manual 
procedure and hence is time consuming and labor intensive. There is currently no 
clinically useful techniques for determining optimal CRT configurations on a 
substantially automatic or a real-time, machine readable basis. 

[0012] Hence, a need arises for a method for automated optimization of multi- 
electrode pacing for cardiac resynchronization. 

SUMMARY 

[0013] One embodiment of the present invention provides a system for 
automatically optimizing CRT procedures using a multi-electrode pacing lead. During 
operation, the system performs a first set of iterations to select one or more satellites on 
one or more pacing leads inserted in a patient. A pacing lead includes a plurality of 
pacing satellites, and a pacing satellite includes a plurality of electrodes that can be 
individually addressed and used for transmitting or detecting electric signals. The system 
then performs a second set of iterations to select one or more electrodes on the selected 
satellites. The system further performs a third set of iterations to select one or more 
timing configurations for pacing signals transmitted through one or more of the selected 
electrodes. 
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[0014] A further embodiment of the present invention provides a graphic user 
interface (GUI) for configuring a cardiac pacing lead that accommodates a plurality of 
pacing satellites. The GUI provides an illustration of the pacing lead, which indicates the 
status of a plurality of electrodes on each pacing satellite. The GUI also provides a 
selection field that allows a user to configure a polarity, rate, amplitude, and pulse width 
of a pacing signal that is to be transmitted to the pacing lead. 

[0015] The present invention has important clinical applications in improving 
results from cardiac ^synchronization therapy (CRT). The dP/dt parameter is a well- 
documented proxy for left ventricular contractility. Use of the present invention allows 
optimization of CRT as evaluated by the dP/dt parameter by about 5-50%, more 
specifically from about 10-20%, and most specifically by about 15%. Similar 
improvements can also be seen using changes in responder claims, cardiac output, and 
synchrony index measures. 

BRIEF DESCRIPTION OF THE FIGURES 
[0016] FIG. 1 illustrates the locations of a number of pacing satellites 
incorporated in a multi-electrode pacing lead, in accordance with an embodiment of the 
present invention. 

[0017] FIG. 2 illustrates an exemplary external view of a number of pacing 
satellites, in accordance with an embodiment of the present invention. 

[0018] FIG. 3 presents a high-level flow chart illustrating the automated 
optimization process using hierarchical iterations, in accordance with an embodiment of 
the present invention. 

[0019] FIG. 4 presents a flow chart illustrating, in conjunction with step 302 of 
FIG. 3, an exemplary implementation of the process of identifying valid satellite 
combinations, in accordance to one embodiment of the present invention. 

[0020] FIG. 5 presents a flow chart illustrating, in conjunction with step 304 of 
FIG. 3, an exemplary implementation of the process of identifying valid electrode 
combinations, in accordance with an embodiment of the present invention. 
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[0021] FIG. 6 presents a flow chart illustrating, in conjunction with step 306 of 
FIG. 3, an exemplary implementation of the process of identifying acceptable timing 
configurations, in accordance with an embodiment of the present invention. 

[0022] FIG. 7 presents a flow chart illustrating, in conjunction with step 602 of 
FIG. 6, an exemplary implementation of the process of performing iterative 
measurements with different timing configurations using a coarse-granularity step length, 
in accordance with an embodiment of the present invention. 

[0023] FIG. 8 presents a flow chart illustrating, in conjunction with step 606 of 
FIG. 6, an exemplary implementation of the process of performing iterative 
measurements with different timing configurations using a fine-granularity step length, in 
accordance with an embodiment of the present invention. 

[0024] FIG. 9 illustrates an electrocardiogram (EGG) signal, measured in a pig's 
heart, in response to a right-ventricle pacing signal and a left-ventricle packing signal, in 
accordance to an embodiment of the present invention. 

[0025] FIG. 10 illustrates measuring the timing difference between two cardiac 
signals using their corresponding relative values, in accordance with an embodiment of 
the present invention. 

[0026] FIG. 1 1 is an exemplary illustration of a three-dimensional plot of 
synchrony indices versus the timing delay between two pacing pulses and the locations of 
the pacing electrodes, in accordance with an embodiment of the present invention. 

[0027] FIG. 12 is an exemplary illustration of a two-dimensional plot of a 
synchrony index versus the timing delay between two pacing pulses, in accordance with 
an embodiment of the present invention. 

[0028] FIG. 13 is an exemplary illustration of a linear display for synchrony 
index, in accordance with an embodiment of the present invention. 

[0029] FIG. 14 is an exemplary illustration of X-Y traces of two sensor signals 
for purposes of displaying synchrony, in accordance with an embodiment of the present 
invention. 

[0030] FIG. 15 illustrates an X-Y trace of left-ventricle pressure versus left- 
ventricle lateral strain based on actual data collected by uncalibrated sensors placed in a 
pig's heart, in accordance with an embodiment of the present invention. 
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[0031] FIG. 16 illustrates an X-Y trace of left-ventricle pressure versus right- 
ventricle freewall strain based on actual data collected by uncalibrated sensors placed in a 
pig's heart, in accordance with an embodiment of the present invention. 

[0032] FIG. 1 7 illustrates an X-Y trace of left- ventricle pressure versus right- 
ventricle septal strain based on actual data collected by uncalibrated sensors placed in a 
pig's heart, in accordance with an embodiment of the present invention. 

[0033] FIG. 1 8 is an exemplary illustration of a radar-scan-like display of two 
sensor signals for indicating synchrony, in accordance with an embodiment of the present 
invention. 

[0034] FIG. 19 is an exemplary illustration of a graphic-equalizer-like display 
of multiple sensor signals for indicating synchrony, in accordance with an embodiment of 
the present invention. 

[0035] FIG. 20 is an exemplary illustration of a graphic-equalizer-like display 
of multiple sensor signals indicating substantial synchrony, in accordance with an 
embodiment of the present invention. 

[0036] FIG. 21 illustrates a measurement of electromechanical delay (EMD) for 
measuring synchrony, in accordance with an embodiment of the present invention. 

[0037] FIG. 22 illustrates a measurement of EMD based on one of multiple 
strain gauge signals for measuring synchrony, in accordance with an embodiment of the 
present invention. 

[0038] FIG. 23 illustrates an exemplary graphic user interface (GUI) for 
configuring a multi-electrode pacing lead and for displaying measured signals, in 
accordance with an embodiment of the present invention. 

[0039] FIG. 24 illustrates another exemplary GUI for configuring a multi- 
electrode pacing lead, in accordance with an embodiment of the, present invention. 

[0040] FIG. 25 is an exemplary illustration of a display showing raw data 
collected from different electrodes in response to different electric fields in an electric 
tomography measurement, in accordance with an embodiment of the present invention. 

[0041] FIG. 26 is an exemplary illustration of a display showing comparison of 
data collected from different electrodes in response to different electric fields in an 
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electric tomography measurement, in accordance with an embodiment of the present 
invention. 

[0042] FIG. 27 illustrates an exemplary GUI for confirming connectivity 
between the pacing system, the pacing leads, and the pacing electrodes, in accordance 
with an embodiment of the present invention. 

[0043] FIG. 28 illustrates an exemplary GUI for logging in and loading saved 
cardiac pacing and measurement data of a patient, in accordance with an embodiment of 
the present invention. 

[0044] FIG. 29 is a Pacing Capture Chart which illustrates an exemplary GUI 
where, Green dots mean Anode, Red dots mean Cathode, and the number in the middle 
of the circle is the satellite number. 

DETAILED DESCRIPTION 

[0045] The following description is presented to enable any person skilled in 
the art to make and use the invention, and is provided in the context of a particular 
application and its requirements. Various modifications to the disclosed embodiments 
will be readily apparent to those skilled in the art, and the general principles defined 
herein may be applied to other embodiments and applications without departing from the 
spirit and scope of the present invention. Thus, the present invention is not intended to 
be limited to the embodiments shown, but is to be accorded the widest scope consistent 
with the principles and features disclosed herein. 

[0046] The data structures and the programs for performing operations 
described in this detailed description are typically stored on a computer readable storage 
medium, which may be any device or medium that can store code and/or data for use by a 
computer system. This includes, but is not limited to, application specific integrated 
circuits (ASICs), field-programmable gate arrays (FPGAs), semiconductor memories, 
magnetic and optical storage devices such as disk drives, magnetic tape, CDs (compact 
discs), and DVDs (digital versatile discs or digital video discs). 

[0047] Embodiments of the present invention provide a system that 
automatically optimizes the configuration of a multi-electrode pacing lead for cardiac 
pacing and for detecting signals from sensors incorporated on the lead. According to one 
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embodiment, a pacing lead includes a number of pacing circuit elements, referred to as 
"satellites," which further include a number of electrodes. Each satellite and the 
electrodes therein can be individually addressed and controlled through a low-power 
control chip. The satellites are coupled to a pacing-control system, which can be a 
5 pacing can or an external control system. One embodiment of such a multi-electrode 
pacing lead is described in a prior application by some of the present inventors. 

[0048] A multi-electrode pacing lead can contain a number of satellites, each of 
which containing multiple electrodes. A large number of combinations are therefore 
possible for transmitting bi-polar or uni-polar (when, for example, the pacing can is used 
1 0 as signal return) pacing signals. Furthermore, identifying the optimal timing delay 

between two pacing signals sent to two locations can involve hundreds or thousands of 
trial data points, because the delay between heartbeats is on the order of one second, and 
the cardiac synchrony is measured on the order of milliseconds. The number of possible 
timing configurations, when multiplied by the number of possible satellite/electrode 
1 5 combinations, can make the cost of optimization prohibitively high. For example, in a 
bi-ventricular pacemaker system with at least one lead having eight satellites where each 
satellite has four electrodes, a clinician can face billions of configurations with different 
satellite/electrode/timing combinations. Using TDI to scan manually through all possible 
configurations is impossible. Without an automated optimization process, the clinician 
20 can only experiment with a small number of configurations. Such an empirical approach 
is not optimally reliable or effective. 

[0049] Embodiments of the present invention provide a system that can quickly 
identify valid satellite and electrode combinations by performing hierarchical iterations. 
The system then uses a directed search to find optimized timing configuration for each 
25 electrode combination. The system subsequently detects peaks of the best synchrony and 
provides the corresponding satellite/electrode and timing configurations to the physician. 
If the identified configurations are unacceptable, the system can further indicate that a 
change of lead placement may be necessary. 

30 Multi-Electrode Pacing Lead 
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[0050] FIG. 1 illustrates the locations of a number of pacing satellites 
incorporated in a multi-electrode pacing lead, in accordance with an embodiment of the 
present invention. A pacing and signal detection system 101 provides extra-cardiac 
communication and control elements for the overall system. In some embodiments, 
pacing and signal detection system 101 may be, for example, a pacing can of a 
pacemaker residing in an external or extra-corporeal location. 

[0051] Right ventricular lead 102 emerges from pacing and signal detection 
system 101 and travels from a subcutaneous location from pacing and signal detection 
system 101 into the patient's body (e.g., preferably, a subclavian venous access), and 
through the superior vena cava into the right atrium. From the right atrium, right 
ventricle lead 102 is threaded through the tricuspid valve to a location along the walls of 
the right ventricle. The distal portion of right ventricular lead 102 is preferably located 
along the intra-ventricular septum, terminating with a fixation in the right ventricular 
apex. Right ventricular lead 102 includes satellites positioned at locations 103 and 104. 
The number of satellites in ventricular lead 102 is not limited, and may be more or less 
than the number of satellites shown in FIG. 1 . 

[0052] Similarly, left ventricular lead 105 emerges from pacing and signal 
detection system 101, following substantially the same route as right ventricular lead 102 
(e.g., through the subclavian venous access and the superior vena cava into the right 
atrium). In the right atrium, left ventricular lead 105 is threaded through the coronary 
sinus around the posterior wall of the heart in a cardiac vein draining into the coronary 
sinus. Left ventricular lead 105 is provided laterally along the walls of the left ventricle, 
which is likely to be an advantageous position for bi-ventricular pacing. FIG. 1 shows 
satellites positioned at locations 106 and 107 along left ventricular lead 105. 

[0053] Right ventricular lead 1 02 may optionally be provided with pressure 
sensor 108 in the right ventricle. A signal multiplexing arrangement allows a lead to 
include such active devices (e.g., pressure sensor 108) for pacing and signal collection 
purposes (e.g., right ventricular lead 102). Pacing and signal detection system 101 
communicates with each of the satellites at locations 103, 104, 106 and 107. The 
electrodes controlled by the satellites may also be used to detect cardiac depolarization 



WO 2006/105474 



PCT/US2006/012246 



11 

signals. Additionally, other types of sensors, such as an accelerometer, strain gauge, 
angle gauge, temperature sensor, can be included in any of the leads. 

[0054] FIG. 2 illustrates an exemplary external view of a number of pacing 
satellites, in accordance with an embodiment of the present invention. According to one 
embodiment, a pacing lead 200 (e.g., right ventricular lead 102 or left ventricular lead 
105 of FIG. 1) accommodates two bus wires SI and S2, which are coupled to a number 
(e.g., eight) of satellites, such as satellite 202. FIG. 2 also shows satellite 202 with a 
zoom-in view. Satellite 202 includes electrodes 212, 214, 216, and 218, located in the 
four quadrants of the cylindrical outer walls of satellite 202. Each satellite also contains 
a control chip which communicates with a pacing and signal-detection system to receive 
configuration signals that determine which of the four electrodes are to be coupled to bus 
wires SI or S2. 

[0055] The configuration signals, the subsequent pacing pulse signals, and the 
analog signals collected by the electrodes can all be communicated through bus wires SI 
and S2, in either direction. Since both digital and analog signals can be communicated 
through bus wires SI and S2, no digital-to-analog or analog-to-digital converters need to 
be provided in the satellites. Although shown in a symmetrical arrangement, electrodes 
212, 214, 216 and 218 may be offset along lead 200 to minimize capacitive coupling 
among these electrodes. The quadrant arrangement of electrodes allows administering 
pacing current via electrodes oriented at a preferred direction, for example, away from 
nerves, or facing an electrode configured to sink the pacing current. Such precise pacing 
allows low-power pacing and minimal tissue damage caused by the pacing signal. 

Hierarchical Iterative Optimization 

[0056] Since a multi-electrode lead can provide a large number of electrode 
combinations for pacing purposes, to identify effective combinations quickly is very 
advantageous. Theoretically, a system can automatically try every possible electrode 
combination and timing configuration in a brute-force manner, which is time consuming. 
Embodiments of the present invention employ three levels of hierarchical iterations to 
expedite the optimization process. With such a process, the system can quickly eliminate 
satellite and electrode combinations that cannot effectively produce a cardiac response. 
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[0057] For the purposes of this application, "valid" means an electrode satellite 
group which captures heart muscle, and in the case of a segmented electrode one which 
also avoids phrenic nerve capture. Also for the purposes of this application, "configure" 
means the timing and electrode selection provided-that is on-off and polarity selection. 

[0058] Because the system aims to identify the valid satellite combinations, 
electrode combinations, and timing configurations, one approach is to conduct the trial 
pacing in three corresponding levels. That is, the system first experiments with all 
possible satellite combinations and identifies valid ones. Valid meaning pacing capture 
occurs and, at the segmented electrode selection level, phrenic nerve capture does not 
occur. Then, based on the valid satellite combinations, the system further experiments 
with all possible electrode combinations and identifies those that are valid. Finally, the 
system experiments with all possible timing configurations and identifies those that are 
optimal. This hierarchical iteration process can also be performed with a specified or 
limited scope of possible configurations. For example, a physician can specify certain 
satellites, electrodes, or timing configurations with which the system is to experiment. In 
general, this hierarchical iteration process can be used to perform both global and local 
optimization. 

[0059] FIG. 3 presents a high-level flow chart illustrating the automated 
optimization process using hierarchical iterations, in accordance with an embodiment of 
the present invention. During operation, the system first identifies the valid satellite 
combinations (step 302). The system then identifies the valid electrode combinations 
among electrodes on the valid satellites identified in step 302 (step 304). Next, the 
system identifies acceptable timing configurations using valid electrodes identified in 
step 304 (step 306). The system subsequently provides the optimal 
satellite/electrode/timing configurations to a physician (step 308). 

[0060] During the first set of iterations, the system experiments with all satellite 
combinations. In one embodiment, the system uses all the electrodes on a satellite when 
pacing through that satellite. When performing the iterations, the system sequentially 
selects satellite pairs, wherein one satellite is used as signal source and another as signal 
sink for bipolar pacing. The system can also include uni-polar pacing in the iterations. 
In uni-polar pacing, the pacing can is used as signal source (or sink) and one satellite is 
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used as signal sink (or source). The uni-polar pacing can be viewed as a special instance 
of bipolar pacing, wherein the can functions as a satellite with only one electrode. In one 
embodiment, the system performs n+1 P 2 iterations, where P denotes the permutation 
operation, n denotes the total number of satellites, and (w+1) indicates that the pacing can 
is considered as a satellite for uni-polar pacing. In a further embodiment, more than one 
satellite can be used as the signal source or signal sink. 

[0061] During each iteration, the system transmits a pulse through the satellite 
pair. If the pulse results in a cardiac response, such as a contraction, the satellite pair is 
considered as a valid combination, which is saved. Otherwise, the system continues to 
experiment with the next untested satellite combination. The presence of valid satellite 
combinations may depend on lead placement. If the system does not identify any valid 
satellite combination, the system can indicate that a change of the lead placement. 

[0062] FIG. 4 presents a flow chart illustrating, in conjunction with step 302 of 
FIG. 3, an exemplary implementation of the process of identifying valid satellite 
combinations, in accordance to one embodiment of the present invention. During 
operation, the system first configures a signal-source satellite (step 402) and a signal-sink 
satellite (step 404). The system then provides a pacing pulse through the configured 
satellites (step 406), and determines whether the pulse results in a cardiac response (step 
408). 

[0063] If a response is present, the system saves the satellite combination as a 
valid combination (step 410). Otherwise, the system proceeds to determine whether all 
signal-sink satellites have been tested in combination with the current signal-source 
satellite (step 412). If not, the system chooses the next untested satellite as the signal- 
sink satellite (step 414) and configures this next satellite (step 404). If all signal-sink 
satellites have been tested, the system further determines whether all signal-source 
satellites have been tested (step 416). If not, the system chooses the next untested 
satellite as the signal-source satellite (step 418) and configures this next satellite (step 
402). Otherwise, the system exits the current set of iterations for selecting valid satellite 
combinations. 

[0064] After identifying valid satellite combinations, the system performs the 
second set of iterations to select valid electrode combinations from the electrodes on all 
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the valid satellites. Thus, by eliminating satellite combinations that do not produce a 
cardiac response, the system also eliminates the associated electrode combinations. This 
elimination significantly reduces the search space for the second set of iterations. 

[0065] FIG. 5 presents a flow chart illustrating, in conjunction with step 304 of 
FIG. 3 5 an exemplary implementation of the process of identifying valid electrode 
combinations, in accordance with an embodiment of the present invention. During 
operation, the system first selects a valid satellite combination (step 501). The system 
then configures a signal-source electrode on the valid signal-source satellite (step 502) 
and a signal-sink electrode on the valid signal-sink satellite (step 504). The system 
subsequently provides a pacing pulse through the configured electrodes (step 506), and 
determines whether the pulse results in a cardiac response (step 508). 

[0066] Some of the inventors' clinical testing indicates that pacing using one or 
two electrodes in a quadrant reduces the capture threshold voltage and increases the 
electrode impedance. This effect is in contrast to when all four quadrants are employed 
in a ring like manner. The capture threshold reduction ranges from about 20% to as 
much as 75%, more specifically from about 30%-60%, and most specifically about 50%. 
The impedance increase is in the range of 100% (1 vs. 4 electrodes) 

[0067] This combination has a powerful effect on power consumed in the 
pacing function given P = i*i*R, lowering V, reducing current. Increasing R lowers 
current. Optimization may focus on finding the one or two optimal electrodes in a 
quadrant. The risk is the optimal electrode may shift until the lead fixates several weeks 
after implant. 

[0068] If a response is present, the system saves the electrode combination as a 
valid combination (step 510). Otherwise, the system proceeds to determine whether all 
electrodes on the signal-sink satellite have been tested in combination with the current 
signal-source electrode (step 512). If not, the system chooses the next untested electrode 
on the signal-sink satellite as the signal-sink electrode (step 514) and configures this next 
electrode (step 504). If all signal-sink electrodes have been tested, the system further 
determines whether all electrodes on the signal-source satellite have been tested (step 
516). If not, the system chooses the next untested electrode on the signal-source satellite 
as the signal-source electrode (step 518) and configures this next electrode (step 502). 
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Otherwise, the system determines whether all valid satellite combinations have been 
tested (step 520). If not, the system proceeds to select the next valid satellite 
combination (step 501). Otherwise, the system exits the current set of iterations for 
selecting valid electrode combinations. 

[0069] After identifying the valid satellite and electrode combinations, the 
system then enters the third set of iterations to identify optimal timing configurations. 
The goal is to find a timing configuration, which can be the time delay between two 
pacing signals or between one pacing signal and the contraction of a specific location, so 
that the heart can attain the best synchrony. Because the synchrony can be sensitive of a 
timing difference of only a few milliseconds while the period between two consecutive 
heartbeats is on the order of one second, hundreds or even over a thousand of test data 
points are present By example, there can be about 100-10,000 test data points, more 
specifically about 500-1,000 test data points, and most specifically about 700 test data 
points One approach to find the timing configuration with the best synchrony is to 
perform a brute-force, exhaustive search by experimenting with every possible timing 
configuration with a fixed step length, for example, 1 ms. This approach, however, can 
be very time consuming and may not be clinically practical. 

[0070] One embodiment of the present invention employs two rounds of 
iteration for searching timing configurations that produce the best synchrony. In the first 
round of iterations, the system scans through all the valid electrode combinations, and, 
for each electrode combination, the system experiments with different timing 
configurations using a coarse-granularity step length, and records the measured 
corresponding synchrony. The term "electrode combination" here refers to a first valid 
electrode configured as a signal source and a second valid electrode configured as a 
signal sink, and hence includes different valid satellite combinations. Based on a 
synchrony- versus-delay plot, the system identifies broad peaks of best synchrony values 
using a given set of criteria. In the second round of iterations, the system focuses on and 
re-scans these broad peaks, and experiments with different timing configurations using a 
fine-granularity step length. The system can thereby quickly and accurately identify 
timing configurations that result in the best available synchrony. 
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[0071] This two-level iterative process can be analogized to the process of 
finding a target object under a microscope. A viewer first uses a low-power objective 
lens to scan the specimen and locates the approximate region where the target is. The 
user then places the target region in the center of the view field and changes to a high- 
power objective lens. By adjusting the focus and fine-tuning the position of the 
specimen, the user can quickly obtain a clear view of the target object. In the present 
system, the first round of iterations with a coarse-granularity step length resembles 
locating the target region with a low-power objective lens; and the second round of 
iterations with a fine-granularity step length resembles the focusing and fine-tvming with 
the high-power objective lens. 

[0072] The ratio between the coarse-granularity and find-granularity step 
lengths can be set at a default value or be determined by a clinician. In one embodiment, 
the coarse-granularity step length is approximately about 5-10 times the fine-granularity 
step length. In a further embodiment, the coarse-granularity step length is about 5-10 ms. 
Most broadly, this step length can be between about 4-40ms, more specifically from 
about 10-3 0ms, and most specifically about 15ms. By contrast, the fine-granularity step 
length is between about 0.05 and 4.00ms, more specifically between about 0.10 and 
2.00ms, and most specifically about 1 ms. Other step-length ratios and values are also 
possible. 

[0073] FIG. 6 presents a flow chart illustrating, in conjunction with step 306 of 
FIG. 3, an exemplary implementation of the' process of identifying acceptable timing 
configurations, in accordance with an embodiment of the present invention. During 
operation, the system performs a first round of iterative measurements of synchrony with 
different timing configurations for all the valid electrode combinations using a coarse- 
granularity step length (step 602). The system subsequently identifies broad peaks for 
synchrony (step 604). 

[0074] Next, the system performs a second round of iterative measurements of 
synchrony with different timing configurations centered around the broad peaks, using a 
fine-granularity step length (step 606). The system then identifies the most desirable 
timing configuration(s) and the corresponding electrode combination(s) (step 608). The 
system subsequently provides a pulse through the optimal electrode combination with the 
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optimal timing configuration (step 610), and measures the synchrony (step 612). The 
system further saves and displays the measured synchrony to a clinician (step 614). 

[0075] FIG. 7 presents a flow chart illustrating, in conjunction with step 602 of 
FIG. 6, an exemplary implementation of the process of performing iterative 
measurements with different timing configurations using a coarse-granularity step length, 
in accordance with an embodiment of the present invention. The system selects a valid 
electrode combination (step 702), and selects an initial pacing delay (step 704). The 
system then provides pacing pulse according to the selected pacing delay through the 
selected electrodes (step 706). Next, the system measures synchrony (step 708), saves 
and displays the measured synchrony (step 710). 

[0076] The system subsequently determines whether all pacing delays have 
been tested (step 712). If not, the system increments the pacing delay by the coarse- 
granularity step length (714) and provides another pacing pulse with the new pacing 
delay (step 706). If all pacing delays have been tested, the system further determines 
whether all valid electrode combinations have been tested (step 716). If not, the system 
selects the next valid electrode combination (step 702) and performs the iterations again. 
Otherwise, the system exits. 

[0077] After the first round of iterations, the system identifies the broad peaks 
of synchrony. In one embodiment, the system plots synchrony versus pacing delay for 
every valid electrode combination, thereby creating a number of two-dimensional 
plottings. On each plotting, the system can first identify the peaks, which are data points 
with the highest synchrony indices. The system then expands the peak to include a 
number of neighboring data points, for example 100 data points, on both sides of the 
peak. These data points define the width a corresponding broad peak, which is 
subsequently used for locating the highest-synchrony. In one embodiment, data points 
immediately adjacent to a data point with the same or a higher synchrony are eliminated, 
since such data points are part of the same broad peak. The system can further trace 
paths for the remaining data points. If the derivative of the path between two data points 
does not show a sign change, then the data points are identified as being part of the same 
broad peak. 
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[0078] In a further embodiment, the width of each broad peak is determined by 
a pre-determined cut-off synchrony value, which can be set at a percentage of the 
synchrony value of the peak point. For example, the cut-off value can be between about 
20% and 90%, more specifically between 35% and 66%, and most specifically about 
50% of the peak synchrony value. 

[0079] In yet a further embodiment, the system compares two identified peaks 
by first multiplying each peak's height with the peak's width at a cut-off threshold, and 
then comparing the corresponding products. This product can be referred to as a cross 
sectional area. The peaks with the largest cross sectional areas or a cross sectional area 
above a threshold is identified as a broad peak. 

[0080] The system can also visually display a synchrony map, which indicates 
the synchrony with different timing and electrode configurations. A clinician can then 
visually identify the synchrony peaks. In one embodiment, the clinician selects a peak 
that is displayed on a computer screen using a mouse, trackball, touch-screen, or joystick. 

[0081] FIG. 8 presents a flow chart illustrating, in conjunction with step 606 of 
FIG. 6, an exemplary implementation of the process of performing iterative 
measurements with different timing configurations using a fine-granularity step length, in 
accordance with an embodiment of the present invention. The system first selects a 
synchrony broad peak (step 802), and increments the pacing delay from the value 
corresponding to the peak by the fine-granularity step length (step 804). The system then 
measures, saves, and displays the synchrony (step 806). The system further determines 
whether all pacing-delay values higher than the value corresponding to the peak have 
been tested (step 808). If not, the system increments the pacing delay by the fine- 
granularity step length (step 809) and continues to measure the new synchrony (step 
806). 

[0082] If all pacing-delay values higher than the value corresponding to the 
peak value have been tested, the system then decrements the pacing delay from the value 
corresponding to the peak by the fine-granularity step length (step 810). The system 
measures, saves, and displays the synchrony (step 812), and determines whether all 
pacing-delay values lower than the value corresponding to the peak have been tested 
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(step 814). If not, the system continues to decrement pacing delay by the fine-granularity 
step length (step 815) and measures the new synchrony (step 812). 

[0083] If all pacing-delay values lower than the value corresponding to the peak 
have been tested, the system further determines whether all the broad peaks have been 
tested (step 816). If not, the system continues to select the next broad peak and repeats 
the iterative process (step 802). Otherwise, the system exits. 

Automated Optimization Using Scripting 

[0084] Typically, the system uses different hardware and middleware modules 
to control the multi-electrode pacing lead. To perform the automated optimization 
process, these modules are driven with different parameters in a given sequence, while 
response data are collected simultaneously. One possible approach to implement such a 
system is to develop customized, embedded software tailored to specific hardware used 
in the system. Such a system, however, can be costly and inflexible, because the control 
software is closely coupled to the underlying hardware and middleware, and cannot be 
easily modified. 

[0085] One embodiment of the present invention uses a scripting language to 
control the hardware and middleware. A scripting language, compared with other types 
of programming languages, are usually simpler to learn and use, and does not need to be 
complied. The scripting language is interpreted by a computer at run-time, so that the 
machine can execute instructions immediately. Using a scripting language allows the 
present system to de-couple the underlying middleware that drives the hardware, from 
the high-level logical controls. For example, the system can include a first middleware 
module that configures the multi-electrode lead and provides pacing pulses with different 
timing configurations, and a second module that collects sensor response signals. The 
system can then perform the automated optimization process using a script program that 
controls these middleware modules to carry out the hierarchical iterations described 
above. 

[0086] Using scripting language has many advantages. Scripting effectively de- 
couples the logical programming from the underlying hardware and middleware. Such 
de-coupling makes the present system modular and allows each hardware or middleware 
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module to be used for more than one application. Scripting also provides a wide range of 
flexibility and allows a user to reconfigure parameters for the optimization easily. With a 
graphic scripting tool, one can create specific scripts by manipulating (e.g., drag and 
drop) graphic objects and by inputting numerical or alphabetical values. Such a tool is 
particularly useful to clinicians who are not familiar with coding computer programs. 
For example, a clinician can use the graphic scripting tool to create a special test 
sequence, which may or may not involve the aforementioned hierarchical iterations, to 
experiment with a specific set of electrode/timing configurations. The graphic scripting 
tool can generate a corresponding script to perform the test sequence. 

[0087] Another advantage of using scripting language is that a user can save and 
load test data with great flexibility. A script can specify the format in which the system 
saves data. The system can also allow a user to load previously saved data obtained for a 
specific patient at a specific time. Furthermore, a script can manage additional complex 
data processing tools to perform various manipulation on the collected data. 

[0088] In one embodiment, the present system includes a configuration and 
driver module and a data collection module, which are based on LabVIEW-compatible 
hardware from National Instruments Inc. of Austin, Texas. 

Synchrony Measurement 

[0089] Embodiments of the present invention provide a number of ways for 
measuring and displaying synchrony information. The various ways of displaying 
synchrony described herein allow a clinician to observe synchrony in real-time and in an 
intuitive manner. Additionally, the input of the clinician can assist with the optimization 
process in, for example, choosing the broad peaks. 

[0090] Typically, the synchrony is measured as the timing difference between 
two sensor signals in response to the pacing signals. FIG. 9 illustrates an 
electrocardiogram (ECG) signal, measured in a pig's heart, in response to a right- 
ventricle pacing signal and a left-ventricle packing signal, in accordance to an 
embodiment of the present invention. The pacing signals include a left-ventricle signal 
902 and a right-ventricle signal 906, both of which provide negative pulses. A right- 
ventricle pulse 908 is followed by a left-ventricle pulse 904, which is further followed by 
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a peak 910 of an ECG signal 912. The timing difference between pulses 904 and 908 
constitutes a specific timing configuration. The system measures the synchrony as a 
function of the pulse-timing differences. 

[0091] Synchrony measurement typically involves measuring the timing 
difference between corresponding characteristic points of two sensor signals, which 
indicate cardiac contractions at two different locations. Embodiments of the present 
invention allow the system to use relative values, instead of absolute values, of sensor 
signals for measuring synchrony. Relative values are based on signals from an 
uncalibrated sensor. The system can thereby avoid calibration of these sensors, which 
can be time consuming and expensive, because the size, location, and signal drifting of 
the sensors can make calibration very difficult. Furthermore, using relative and 
dimensionless values reduces the processing required for conversion or normalization. 
Consequently, the response time for synchrony measurement can be reduced, which in 
turn expedites the optimization process. In further embodiments, the system may 
normalize or weight the sensor signals before measuring the synchrony. 

[0092] Furthermore, embodiments of the present invention can tolerate signal 
drift without using a filter, because signal drift does not significantly change the position 
of a peak. Nevertheless, an optional low frequency filter can be used to remove drift. For 
example, the hear beats at a frequency of approximately 1 Hz. The measured sensor 
signals therefore also have a frequency of about 1 Hz. A typical drift is on the order of 
0.1 to 0.2 Hz. A high-pass filter may have a cut-off frequency from about 0.05 to 0.2 Hz, 
more specifically between 0.10 and 0.15 Hz, and most specifically at about 0.1 Hz can be 
used to remove the drift. 

[0093] FIG. 10 illustrates measuring the timing difference between two cardiac 
signals using their corresponding relative values, in accordance with an embodiment of 
the present invention. A strain gauge in the left ventricle generates a signal 1002, and a 
second sensor in the right ventricle generates a signal 1004. The characteristic points 
chosen for comparison are the maximum-slew-rate points of the two signals. Slew rate is 
the rate of change (or slope) of a signal. The maximum slew rate point of a signal is 
where the signal changes at the highest rate. 
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[0094] The system determines the synchrony by measuring the timing 
difference, At, between the maximum-slew-rate points of the two signals. In one 
embodiment, the system calculates a synchrony index based on the absolute value of At. 
A better synchrony is attained with a smaller absolute value of At. The system can 
further group the measured synchrony into three groups, namely "good," "OK," and 
"bad," based on two thresholds x and y. If At<x, then the synchrony is designated as 
"good." If x<At < y, the synchrony is designated as "OK." If At>y, the synchrony is 
designated as "bad." 

[0095] In a further embodiment, the system can calculate \At T - At\, the absolute 
value of the difference between At and a target time delay At T . The optimal synchrony 
corresponds to the smallest absolute value of the difference. In a preferred embodiment, 
the target time delay is zero. 

[0096] Further embodiments of the present invention can employ time-average 
operations to smoothen the sensor signals and to eliminate interferences due to noises 
and distortions. The system can automatically configure a time window for averaging the 
signals based on the instant rate of change of the signal. If the instant rate of change is 
high and the sign of the signal slope varies frequently, the signal may experience heavy 
noise and interferences, and the system can use a larger averaging time window. The 
system can also allow a user to specify the averaging time window. Although performing 
time averaging to the signal precludes a real-time measurement of the synchrony, doing 
so can make the measurements more accurate. 

[0097] In a further embodiment, the system computes a cross-correlation 
between the two sensor signals. Computing the cross-correlation allows use of the entire 
signal, instead of just one point, to measure synchrony. Doing so helps reduce errors 
caused by noise, interferences, or distortions, which may be due to harmonics. For 
example, the system can compute a cross-correlation function based on left- ventricle 
signal 1002 and right-ventricle signal 1004. The cross-correlation of two sampled (time- 
discrete) signals can be computed as: 

]T x-x • y-y 

a/E(^^) 2 J-E(^^) 2 J 

wherein x and 3^ are the amplitude of the two signals, respectively, and wherein -l<r<l is 
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the cross-correlation coefficient. When the cross-correlation coefficient r=l, a complete 
positive correlation is present between the two signals. That is, the signals are in phase, 
and the best synchrony is attained. When r=-l, there is complete negative correlation 
between the signals, and the signals are completely out of phase. The value of r indicates 
the degree to which the correlation lies between full synchrony and complete negative 
synchrony. When r=0 5 the signals are considered uncorrelated. 

[0098] In one embodiment, the system applies thresholds, x and y, to determine 
the quality of the synchrony. If r>x, then the synchrony is designated as "good." If 
x>r>y, the synchrony is designated as "OK." If y>r, the synchrony is designated as 
"bad." In a further embodiment, thresholds may designate results in a range bounded on 
both sides of zero as "good," values adjacent to either side of the range as "OK," and 
values further away from either side of the rage as "bad." 

[0099] Because the system allows use of relative, or dimensionless, values 
signals from different types of sensors can be used to measure synchrony. For example, 
the two signals shown in FIG. 10 can be from a strain gauge and an accelerometer, 
respectively. In general, any sensors capable of indicating cardiac contraction can be 
used for synchrony measurements. Such sensors include, but are not limited to, angle 
gauges, strain gauges, pressure sensors, electropotential sensors, volume flow sensors, 
piezoelectric sensors, and accelerometers. Furthermore, a sensor can be incorporated on 
a signal chip which is included in the lead. 

Synchrony Display 

[00100] Embodiments of the present invention employ a number of ways to 
display or indicate synchrony information, so that a clinician can observe the response 
signals and assist with the optimization process. 

[00101] FIG. 1 1 is an exemplary illustration of a three-dimensional plot of 
synchrony indices versus the timing delay between two pacing pulses and the locations of 
the pacing electrodes, in accordance with an embodiment of the present invention. The 
X -axis of the display indicates the timing delay between two pacing signals, one of which 
is transmitted to the left ventricle, and the other to the right ventricle. In this example, 
one of the pacing electrodes is fixed, and the location of the other one is varied by 
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selecting different electrodes on a multi-electrode lead. The Y-axis indicates the location 
of the varying electrode. The Z-axis indicates a synchrony index, which can be a cross- 
correlation coefficient between signals from two sensors placed in the two ventricles. 
The display illustrates two three-dimensional peaks, 1 102 and 1 104. Peak 1 102 is lower 
5 than peak 1 1 04. The system may select peak 1 1 04 as the broad peak, because peak 1 1 02 
is too low or too narrow. 

[00102] FIG. 12 is an exemplary illustration of a two-dimensional plot of a 
synchrony index versus the timing delay between two pacing pulses, in accordance with 
an embodiment of the present invention. The X-axis indicates the timing delay between 

10 two pacing signals, one of which is transmitted to the left ventricle, and the other to the 
right ventricle. The Y-axis indicates a synchrony index. Note that this two-dimensional 
display only plots different synchrony-index values at different timing configurations 
with a given electrode combination. Such a two-dimensional display can be considered 
as one of the possible cross-sections of a three-dimensional display as illustrated in FIG. 

15 11. The two-dimensional display illustrates two peaks, 1202 and 1204. Peak 1204 is 

lower than peak 1202. The system may select peak 1202 as the broad peak, because peak 
1204 is too low or too narrow. 

[00103] FIG. 13 is an exemplary illustration of a linear display for synchrony 
index, in accordance with an embodiment of the present invention. This display 

20 indicates the timing delay between the corresponding points of two signals by changing 
the position of an indicator 1308 along a vertical bar. The vertical bar includes three 
regions, 1302, 1304, and 1306. Region 1306 is for timing delays between 0 and 60 ms, 
region 1304 is for delays between 60 and 120 ms, and region 1302 is for delays longer 
than 120 ms. The numerical ranges, such as 60 ms and 120 ms, are for illustrative 

25 purposes only. Other ranges can also be used. In one embodiment, regions 1302, 1304, 
and 1306 can be color-coded. For example, region 1206 can be in green to indicate 
"good" synchrony, region 1304 can be in yellow to indicate "OK" synchrony, and region 
1302 can be in red indicating "bad" synchrony. Indicator arrow 1308 can move in real 
time and indicate the value of the current measured synchrony. A readout value 1310 

30 indicates the numerical value of the timing delay. 
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[00104] FIG. 14 is an exemplary illustration of X-Y traces of relative values of 
two strain sensor signals for purposes of displaying synchrony, in accordance with an 
embodiment of the present invention. On the upper half of FIG. 14 is a real-time X-Y 
plot of two sensor signals. The X-axis indicates a measured strain in the left ventricle, 
and the Y-axis indicates a measured strain in the right ventricle. If the contraction in the 
left ventricle is in complete synchrony with the contraction in the right ventricle, the 
resulting trace would be approximately a line 1404 angled in the upper-right direction, 
because the minimum and the maximum of both signals are reached at approximately the 
same time. If the two signals are out of synchrony to a certain degree, the resulting trace 
becomes more circle or oval like, such as trace 1402. If the two signals are completely 
out of synchrony and demonstrate a 180-degree phase difference, the resulting trace 
becomes approximately a line angled in the upper-left direction (not shown here). 

[00105] Generally, the measured signals can include noise and distortion, and the 
resulting trace appears to be somewhere between a line and an oval. The more the trace 
resembles a line in the upper right direction, the better the synchrony. This display 
provides a clinician with a real-time, intuitive visual display of the synchrony without 
any quantitative processing. Such a qualitative presentation of synchrony does not 
require the use of absolute or calibrated values of sensor signals. The system uses 
relative and dimensionless values of two strain sensors to produce the trace. In addition, 
the system uses the entire signal, instead of just one point, to indicate synchrony, which 
helps reduce error caused by noise in the signal. 

[00106] The X-Y real-time trace display can also display two calibrated signals, 
or a calibrated signals versus an uncalibrated signal. The lower half of FIG. 14 shows an 
X-Y trace of a calibrated left-ventricle pressure signal and an uncalibrated, relative right- 
ventricle strain signal. If the contraction in the left ventricle is in complete synchrony 
with the contraction in the right ventricle, the resulting trace would be approximately a 
line 1408 angled in the upper-right direction. If the two signals are out of synchrony, the 
resulting trace becomes more circle or oval like, such as trace 1406. In this example, the 
pressure measurement is absolute, calibrated, and dimensional (in mmHg), but the strain 
measurement is relative and dimensionless. 
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[00107] FIG. 15-17 illustrate actual real-time X-Y traces with various sensor 
signals measured in a pig's heart. FIG. 15 is a trace of the left-ventricle pressure versus 
left-ventricle lateral strain, FIG. 16 is a trace of the left-ventricle pressure versus right- 
ventricle freewall strain, and FIG. 17 is a trace of the left-ventricle pressure versus the 
right-ventricle septal strain. 

[00108] FIG. 18 is an exemplary illustration of a radar-scan-like display of two 
sensor signals for indicating synchrony, in accordance with an embodiment of the present 
invention. This display provides a sweeping arm 1 808 that rotates one revolution for 
each heartbeat cycle. The display includes an origin, which is the intersection of two 
orthogonal axis's, 1802 and 1804. Sweeping arm 1808 rotates around the origin. As 
sweeping arm 1808 sweeps in real time, signal amplitude is indicated by a radial distance 
from the origin. The instant signal amplitude is traced along sweeping arm 1808 with a 
corresponding radial distance. Magnitude guides 1 806 are provided to allow a user to 
more easily comprehend the magnitude of the signals. 

[00109] Sweeping arm 1808 can trace multiple signals simultaneously. The 
exemplary display in FIG. 18 shows a first trace pattern 1810 from a first signal and a 
second trace pattern 1812 from a second signal. Patter 1810 and pattern 1812 overlap, 
forming an overlapping region 1814, which is shaded. The area of region 1814 indicates 
the degree to which the two signals move in synchrony. In one embodiment, trace 
patterns 1810 and 1812 are each color coded, and overlapping region 1814 is highlighted 
to facilitate easy observation. In a further embodiment, the signals being traced are first 
normalized for accurate comparison. Normalization, however, can be optional so that 
less processing is required. 

[00110] FIG. 19 is an exemplary illustration of a graphic-equalizer-like display 
of multiple sensor signals for indicating synchrony, in accordance with an embodiment of 
the present invention. A horizontal axis 1910 represents a reference time, which in one 
embodiment can be the beginning time of the QRS complex. QRS complex refers to the 
principal deflection in the ECG, representing ventricular depolarization. This equalizer- 
like display includes four channel indicators, 1902, 1904, 1906, and 1908. Each channel 
represents a sensor signal. For example, each channel can represent a signal from one of 
four strain sensors. Further, one channel can represent a strain sensor signal, a second 
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channel can represent a pressure sensor signal, a third channel can represent a volume 
flow meter signal, and a fourth channel can represent an accelerometer signal. These 
channel indicators resemble a graphic equalizer indicator, which is commonly present on 
consumer audio equipment, making the display intuitive to an observer. 

[00111] The bottom of a channel indicator corresponds to a first characteristic 
event in the signal represented by the channel indicator. An example of such first 
characteristic event is the signal maximum. In one embodiment, the first characteristic 
event occurs at approximately the same time as the beginning of the QRS complex. The 
top of each channel indicator represents a second characteristic event in the signal. An 
example of such second characteristic event is the signal minimum. When the heart is in 
good synchrony, different channel indicators, which represent different sensor signals, 
are substantially even with each other. The four signals represented in FIG. 19 are not in 
good synchrony, because the bottom of these indicators are not close to the beginning of 
the QRS complex, and their top are not even with each other. In contrast, the four signals 
shown in FIG. 20 demonstrate good synchrony, because the four channel indicators are 
substantially even with each other. 

[00112] The graphic-equalizer-like display can also operate in conjunction with 
an audio device. The audio device uses atone to indicate the degree of synchrony, 
wherein a higher frequency of the tone indicates better synchrony. In a further 
embodiment, the audio device can create a beating sound for each channel to indicate the 
degree synchrony. 

Measuring Electromechanical Delay (EMD) 

[00113] Embodiments of the present invention measure the electromechanical 
delay (EMD) between an ECG signal and a sensor signal to evaluate the cardiac 
synchrony. In general, a finite period exists between the start of an electrical activity and 
the start of a mechanical event in response. This period is referred to as EMD. EMD can 
be used to indicate the degree of synchrony in a heart. 

[00114] FIG. 21 illustrates a measurement of EMD for measuring synchrony, in 
accordance with an embodiment of the present invention. In this example, a left- 
ventricle distal strain sensor, such as sensor 107 as is shown in FIG. 1, generates a signal 
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2102. A second more proximal strain sensor in the left ventricle, such as the third strain 
sensor 106 as is shown in FIG. 1, generates a second signal 2104. Both signals 2102 and 
2104 are uncalibrated. The EMD is defined as a time delay between a characteristic 
point of an ECG signal 2106 and a characteristic point on one of the sensor signals. In 
this case, the characteristic point of ECG signal 2106 is the start of the QRS complex. 
The characteristic point of a signal is the periodical signal peak. The time difference 
between the two characteristic points, At, is the EMD for the signal. 

[00115] The system can measure EMD by subtracting the starting time of the 
QRS complex from the peak time of a sensor signal. The system can further apply 
thresholds to determine quality of the synchrony. For example, the system may use a set 
of thresholds x and y. For At<x, the synchrony is designated as "good;" for x<A?<y, the 
synchrony is designated as "OK;" and for y<At, the synchrony is designated as "bad." 

[00116] FIG. 22 illustrates a measurement of EMD based on an uncalibrated 
strain gauge signal from a pig's heart, in accordance with an embodiment of the present 
invention. 

Graphic User Interface 

[00117] FIG. 23 illustrates an exemplary graphic user interface (GUI) for 
configuring a multi-electrode pacing lead and for displaying measured signals, in 
accordance with an embodiment of the present invention. This GUI can be used for 
configuring electrodes and timing pace signals transmitted over a multiple electrode lead 
2354. A user can interact with a pacing control system through this GUI using a 
keyboard, mouse, or touch screen. This GUI also allows the user to select and observe 
different dyssynchrony and performance indices to be optimized. 

[00118] The dyssynchrony indices can include intraventricular and 
interventricular dyssynchrony indices. The intraventricular dyssynchrony indices can 
further include a Proteus Time to Systolic Contraction (PTs) index, which is the time 
interval between onset of QRS complex of an ECG signal and onset of systolic initiation. 
Typically, a low PTs value indicates good synchrony. A numeric display 2328, which 
shows an exemplary value of 43 ms, is provided for this time interval. The 
intraventricular dyssynchrony indices can indicate the time of contraction of the left 
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ventricle relative to the septum, which can be measured by using leads 105 and 102 as is 
shown in FIG. 1 . The system can use various sensors, including strain gauges, to obtain 
this measurement. 

[00119] A second intraventricular dyssynchrony index is the Proteus Time to 
Systolic contraction Standard (PTs-SD). After measuring the PTs indices of a number of 
locations, the system can compute the standard deviation of these measurement as an 
index for dyssynchrony. Low PTs-SD values indicate good synchrony. The GUI also 
provides a numeric display 2318 which shows the numeric value of this index. 

[00120] The GUI further provides an interventricular dyssynchrony index, which 
presents the timing delay between a global indicator for the right ventricle and a global 
indicator for the left ventricle (RV-LV). The GUI provides a numeric display 2320 to 
show the numeric value thereof. A global indicator can be the signal from a pressure 
sensor or other types of sensors. 

[00121] The GUI also provides performance indices. The first performance 
index is the blood tissue ratio. To measure the blood tissue ratio, the system can provide 
two electric fields, one at a low frequency and the other at a high frequency. Because 
blood and tissue have different impedance depending on the frequency of the electric 
field, the system can determine the relative amount of blood and tissue between two 
points. When the heart contracts, blood is ejected, the amount of blood decreases, and 
the tissue thickens. Consequently, the blood to tissue ratio decreases. The GUI provides 
a numeric display 2322 to show the numeric value of the blood to tissue ratio. 

[00122] A second performance index is the end systolic volume. In one 
embodiment, the system uses electrodes to both generate and detect electric fields to 
measure approximate systolic volume in real time. A numeric display 2324 provides the 
numeric value of the systolic volume. 

[00123] The GUI also provides an EGG index. A numeric display 2326 shows 
the numeric value of the QRS width. 

[00124] In the illustrated example, for all the indices shown on the GUI, the 
lower the index value the better the performance of the heart. Therefore, the goal is to 
minimize these indices. 
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[00125] The GUI provides a three-dimensional illustration of a multi-electrode 
lead 2354. In this example, lead 2354 includes eight satellites, and each satellite includes 
four electrodes, such as El, E2, E3, and E4 on the first satellite. Other numbers of 
satellites and electrodes are also possible. The GUI further provides an automatic 
optimization button 2327, labeled "GO." When GO button 1748 is selected, the system 
paces through each satellite with all four electrodes. Pulses transmitted through each 
satellite provide sufficiently high voltage for cardiac pacing. 

[00126] The GUI uses a bar graph to display the selected dyssynchrony and 
performance indices for a baseline measurement (without pacing) and for when each 
satellite is pulsed. In one embodiment, the bar graphs are aligned with their 
corresponding satellites. In the left-most bar graph, which display baseline 
measurements, a first bar 2328 indicates the PTs value, a second bar 2330 indicates the 
PTs-SD value, and a third bar 2332 indicates the blood tissue ratio. Similarly, in a 
second graph for measuring indices when a pacing signal is applied to satellite 1, a first 
bar 2334 indicates the PTs value, a second bar 2336 indicates the PTs-SD value, and a 
third bar 2338 indicates a blood tissue ratio. Bar graphs are generated for all eight 
satellites. In this example, satellite 3 exhibits the lowest selected indices, which indicates 
the least dyssynchrony and therefore is the optimal pacing location. 

[00127] A physician may notice the presence of phrenic capture during the 
pacing of satellite 3. In such a case, the physician may select a phrenic capture indicator 
2352. In response, the system cycles through all the electrodes or combination of 
electrodes on satellite 3 to find an electrode or a combination of electrodes that does not 
cause phrenic nerve capture while still pacing the heart. For example, the four electrodes 
on satellite 3 are able to provide at least eight combinations that define eight 
circumferential directions for pacing. These eight electrode combinations may be 
individual electrodes paired with adjacent electrodes, which can be used to obtain a field 
in the direction between the electrode pair. 

[00128] The GUI also provides displays for various selected signals. For 
example, a first graph 2304 is provided for baseline measurements and a second graph 
2314 is provided for measurements corresponding to satellite 3. In graph 2304, a first 
plot 2312 indicates the septal wall motion, a second plot 2312 indicates the left-ventricle 
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wall motion, a third plot 2308 indicates the right-ventricle pressure, and a fourth plot 
2306 indicates the ECG. Pull-down menus 2302 are provided to change the graphs to 
display dyssynchrony from other satellites and to allow display of other indices. 

[00129] FIG. 24 illustrates another exemplary GUI for configuring a multi- 
electrode pacing lead, in accordance with an embodiment of the present invention. This 
GUI provides a two-dimensional illustration of a multi-electrode pacing lead 2402. Also 
illustrated on lead 2402 are the satellites and the four electrode, denoted as E1-E4, on 
each satellite. 

[00130] A user can select a particular electrode by clicking on that electrode. In 
one embodiment, the state of the electrode changes with each click. For example, with 
one click, the electrode is to be coupled to a cathode; with two clicks, the electrode is to 
be coupled to the anode; and with three clicks, the electrode is to be unconnected. The 
state of the electrode can be indicated by its color, and a color legend 2412 indicates the 
correlation between a color and a connection state. 

[00131] The GUI also provides four selection fields for configuring pacing 
parameters. A first selection field 2404 allows the user to specify polarity of the pacing 
signal. A second selection field 2406 allows the user to specify the pacing rate (bpm). A 
third selection field 2408 allows the user to specify the amplitude of the pacing signal. A 
fourth selection field 2410 allows the user to specify the pulse width of the pacing signal. 

[00132] A switch button 2414, when clicked, changes the connection state of the 
pacing lead based on the electrode selections made by the user. A pace button 2416, 
when clicked, commences the pacing based on the parameters provided by the user. A 
clear button 2418, when clicked, stops the pacing and resets the connection state. 

[00133] The GUI also provides a region 2020 for displaying intra-cardiac 
electrogram (iEGM). In addition, the GUI provides a measure button 2022, which, when 
clicked, measures certain parameters of the pacing system, such as the battery voltage, 
the lead impedance, and the lead temperature. Also provided with the GUI is a "Run 
Script File" button 2024. When this button is clicked, the system allows the user to 
specify a script file, based on which the system automatically executes a series of pacing 
and testing with parameters given in the script file. 
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[00134] Embodiments of the present invention also provide a GUI for displaying 
electric signals detected through different electrodes in an electric tomography 
application. In an electric tomography application, an electric field is applied to the 
heart. This field can be provided by, for example, two electrodes on a lead or by an 
external source. Other electrodes on the lead can then pick up a signal induced by the 
field, wherein the amplitude of the signal depends on the location of the detecting 
electrode. Therefore, the induced signal on an electrode indicates the location and 
movement of the electrode. 

[00135] FIG. 25 is an exemplary illustration of a display showing raw data 
collected from different electrodes in response to different electric fields in an electric 
tomography measurement, in accordance with an embodiment of the present invention. 
The GUI as illustrated in FIG. 25 displays raw data of field-induced signals detected by 
electrodes placed at various locations. The display also includes different signals 
induced by different fields. 

[00136] FIG. 26 is an exemplary illustration of a display showing comparison of 
data collected from different electrodes in response to different electric fields in an 
electric tomography measurement, in accordance with an embodiment of the present 
invention. The display also allows the user to select different signals, such as the EKG 
and iEGM, to be displayed simultaneously for comparison. 

[00137] FIG. 27 illustrates an exemplary GUI for confirming connectivity 
between the pacing system, the pacing leads, and the pacing electrodes, in accordance 
with an embodiment of the present invention. The GUI provides an illustration of the 
wires or cables coupling the pacing leads to the pacing-control system, which, in this 
example, is an external computer. The GUI also provides a "check connection" button 
2702, which, when clicked, issues a command to check the illustrated connectivity. If a 
connection appears to have failed, the GUI can issue a warning. In one embodiment, this 
warning is represented by a change of color of the failed wire or cable. 

[00138] FIG. 28 illustrates an exemplary GUI for logging in and loading saved 
cardiac pacing and measurement data of a patient, in accordance with an embodiment of 
the present invention. A physician can specify a patient's identifier, as well as the time 
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and date when the data was stored. In a further embodiment, the log-in GUI can also 
provide a password field for increased security. 

[00139] Embodiments of the present invention can be used in various systems. 
Such systems may include various types of sensors. Such sensors and systems have been 
described in various applications by some of the present inventors. These applications 
also describe multiplexing systems previously developed by some of the present 
inventors with which the present invention can be employed. These applications include: 
U.S. Patent Application No. 10/734490 published as 20040193021 titled: " Method And 
System For Monitoring And Treating Hemodynamic Parameters"; U.S. Patent 
Application No. 1 1/219,305 published as 20060058588 titled: "Methods And Apparatus 
For Tissue Activation And Monitoring"; International Application No. 
PCT/US2005/046815 titled: " Implantable Addressable Segmented Electrodes"; U.S. 
Patent Application No. 1 1/324,196 titled " Implantable Accelerometer-Based Cardiac 
Wall Position Detector"; U.S. Patent Application No. 10/764,429, entitled "Method and 
Apparatus for Enhancing Cardiac Pacing," U.S. Patent Application No. 10/764,127, 
entitled "Methods and Systems for Measuring Cardiac Parameters," U.S. Patent 
Application No. 10/764, 125, entitled "Method and System for Remote Hemodynamic 
Monitoring"; International Application No. PCT/ US2005/0468 1 5 titled: " Implantable 
Hermetically Sealed Structures"; U.S. Application No. 1 1/368,259 titled: " Fiberoptic 
Tissue Motion Sensor"; International Application No. PCT/US2004/041430 titled: " 
Implantable Pressure Sensors,"; U.S. Patent Application No. 1 1/249,152 entitled 
"Implantable Doppler Tomography System,"and claiming priority to: U.S. Provisional 
Patent Application No. 60/617,618; International Application Serial No. 
PCT/USUS05/39535 titled " Cardiac Motion Characterization by Strain Gauge". These 
applications are incorporated in their entirety by reference herein. 

[00140] The foregoing descriptions of embodiments of the present invention 
have been presented for purposes of illustration and description only. They are not 
intended to be exhaustive or to limit the present invention to the forms disclosed. 
Accordingly, many modifications and variations will be apparent to practitioners skilled 
in the art. Additionally, the above disclosure is not intended to limit the present 
invention. The scope of the present invention is defined by the appended claims. 
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What Is Claimed Is: 

1 . A method for performing automatic optimization of cardiac 
^synchronization therapy (CRT), the method comprising: 

5 performing a first set of iterations to select one or more satellites on one or more 

pacing leads inserted in a patient, wherein a pacing lead includes a plurality of pacing 
satellites, and wherein a pacing satellite includes a plurality of electrodes that can be 
individually addressed and used for transmitting or detecting electric signals; 

performing a second set of iterations to select one or more electrodes on the 
10 selected satellites; and 

performing a third set of iterations to select one or more timing configurations for 
pacing signals transmitted through one or more of the selected electrodes. 

2. The method of claim 1, wherein the pacing lead further includes one or 
more sensors. 

15 3. The method of claim 2, wherein the sensor can be an accelerometer, a 

strain gauge, an angle gauge, a pressure sensor, a temperature sensor, or an electrode 
configured to detect a signal induced by an electric field. 

4. The method of claim 1 , wherein performing the first set of iterations 
involves: 

20 iteratively transmitting a pacing signal through the pacing satellites, wherein each 

iteration involves at least one different satellite from the previous iteration; and 

selecting the one or more satellites based on a cardiac response resulting from the 
pacing signal transmitted through such satellites. 

5. The method of claim 4, wherein transmitting the pacing signal through the 
25 satellites involves using at least one satellite or a can as a signal source, and using at least 

one satellite or the can as a signal sink. 
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6. The method of claim 4, wherein transmitting the pacing signal through the 
satellite involves using all the electrodes on the satellite for transmission. 

7. The method of claim 4, wherein selecting the one or more satellites based 
on the cardiac response involves determining whether the pacing signal results in a 
cardiac contraction. 1 

8. The method of claim 1 ? wherein performing the second set of iterations 
involves: 

iteratively transmitting a pacing signal through the electrodes on the selected 
satellites, wherein each iteration involves a different electrode from the previous 
iteration; and 

selecting the one or more electrodes based on a cardiac response resulting from 
the pacing signal transmitted through such electrodes. 

9. The method of claim 8, wherein transmitting the pacing signal through the 
electrodes involves using at least one electrode on one of the selected satellites or a can 
as a signal source, and using at least one electrode on one of the selected satellites or the 
can as a signal sink. 

10. The method of claim 8, wherein selecting one or more electrodes based on 
the cardiac response involves selecting the electrodes through which the pacing signal 
resulting in a cardiac contraction is transmitted. 

1 1 . The method of claim 1 , wherein performing the third set of iteration 
involves: 

iteratively transmitting pacing signals through at least one selected electrode, 
wherein each iteration involves a different timing configuration of the pacing signal from 
the previous iteration; and 

selecting one or more pacing-signal timing configurations based on a cardiac 
response resulting from the pacing signal. 
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1 2. The method of claim 1 1 , 

wherein transmitting the pacing signal through the selected electrodes involves: 

transmitting a first pulse through a first selected electrode at a first 
location, and 

transmitting a second pulse through a second selected electrode at a 
second location; and 

wherein the timing configuration of the pacing signal includes a time delay 
between the first pulse and the second pulse. 

1 3 . The method of claim 1 1 9 wherein iteratively transmitting pacing signals 
through the selected electrodes involves performing a set of coarse-granularity iterations 
and a set of fine-granularity iterations; and 

wherein the timing-configuration difference between the coarse-granularity 
iterations is larger than that between the fine-granularity iterations. 

14. The method of claim 13, wherein the timing-configuration difference 
between the coarse-granularity iterations is approximately 5 to 10 times of that between 
the fine-granularity iterations. 

15. The method of claim 14, wherein the timing-configuration difference 
between the coarse-granularity iterations is at least twice that between the fine- 
granularity iterations. 

16. The method of claim 15, wherein the timing-configuration difference 
between the coarse-granularity iterations is approximately 5 to 10 times of that between 
the fine-granularity iterations. 

1 7. The method of claim 1 6, wherein the timing-configuration difference 
between the coarse-granularity iterations is approximately 5 to 10 ms; and 

wherein the timing-configuration difference between the fine-granularity 
iterations is approximately 1 ms. 
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18. The method of claim 13, wherein selecting the pacing-signal timing 
configurations based on the cardiac response involves measuring a synchrony 
corresponding to the pacing-signal timing configuration; and 

wherein after the coarse-granularity iterations, selecting the pacing-signal timing 
configurations further involves determining a collection of timing configurations that 
produce better synchrony than other timing configurations; and 

wherein the fine-granularity iterations are performed based on the collection of 
timing configurations determined after the coarse-granularity iterations. 

1 9. The method of claim 1 8, wherein determining the collection of timing 
configurations involves identifying a broad peak based on a synchrony-versus-timing 
plot. 

20. The method of claim 1 8 5 wherein determining the collection of timing 
configurations involves grouping the timing configurations based on pre-determined 
threshold. 

2 1 . The method of claim 1 8, wherein measuring the synchrony involves 
measuring a time difference between two response signals indicating cardiac contraction 
at two corresponding locations. 

22. The method of claim 21 , wherein measuring the time difference between 
the two response signals involves comparing corresponding highest-slew-rate points of 
the two signals, thereby allowing relative values of the two signals to be used. 

23. The method of claim 1 8, wherein measuring the synchrony involves 
computing a cross-correlation of two response signals indicating cardiac contraction at 
two corresponding locations. 

24. The method of claim 1 8, wherein measuring the synchrony involves 
comparing electromechanical delays of two response signals indicating cardiac 
contraction at two corresponding locations. 
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25. The method of claim 24, wherein the electromechanical delay of a signal 
is a time delay between a characteristic point of an electrocardiogram (ECG) and a 
characteristic point of the signal. 

26. The method of claim 25, wherein the characteristic point of the ECG is 
the start of the QRS complex; and 

wherein the characteristic point of the signal is the corresponding peak of the 

signal. 

27. The method of claim 1 8, wherein selecting the pacing-signal timing 
configurations involves displaying an illustration of the synchrony corresponding to the 
timing configuration which is viewable to a physician. 

28. The method of claim 27, 

wherein displaying the illustration of the synchrony involves, tracing a point in 
real time on an X-Y coordinate plane; 

wherein the X coordinate of the point indicates a value corresponding to the 
amplitude of a first response signal indicating cardiac contraction at a first location; and 

wherein the Y coordinate of the point indicates a value corresponding to the 
amplitude of a second response signal indicating cardiac contraction at a second location. 

29. The method of claim 27, wherein displaying the illustration of the 
synchrony involves: 

tracing a first point on a polar coordinate plane, wherein the radial coordinate of 
the first point indicates a value corresponding to the amplitude of a first signal indicating 
cardiac contraction at a first location, and wherein the angular coordinate of the first 
point is configured to change at a constant rate; 

tracing a second point on the same polar coordinate plane, wherein the radial 
coordinate of the second point indicates a value corresponding to the amplitude of the 
second signal indicating cardiac contraction at a second location, and wherein the angular 
coordinate of the second point is configured to change at the same constant rate; and 
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highlighting the area enclosed by both traces of the first point and the second 

point. 

30. The method of claim 1, further comprising using scripting language to 
automatically perform the three sets of iterations. 

3 1 . The method of claim 30, further comprising saving values for response 
signals during each set of iteration using scripting language. 

32. The method of claim 30, further comprising loading the saved values 
using scripting language. 

33. The method of claim 30, further comprising saving satellite, electrode, 
and timing configurations using scripting language. 

34. The method of claim 33, further comprising loading the saved satellite, 
electrode, and timing configurations using scripting language. 

35. The method of claim 30, further comprising allowing a physician to create 
pacing procedures by specifying satellite, electrode, and timing configurations. 

36. The method of claim 35, wherein allowing the physician to create pacing 
procedures involves generating scripts that automatically perform the pacing procedures. 

37. A graphic user interface (GUI) for configuring a cardiac pacing lead that 
accommodates a plurality of pacing satellites, the GUI comprising: 

an illustration of the pacing lead, which indicates the status of a plurality of 
electrodes on each pacing satellite; and 

a selection field that allows a user to configure a polarity, rate, amplitude, and 
pulse width of a pacing signal that is to be transmitted to the pacing lead. 
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38. The GUI of claim 37, wherein the illustration of the pacing lead is 
configured to allow the user to configure a connection state and a polarity of each 
electrode. 

39. The GUI of claim 38, wherein the illustration of the pacing lead is 
configured to allow the user to toggle the status of each electrode by clicking on the 
image of the electrode. 

40. The GUI of claim 37, further comprising a field that indicates a battery 
voltage for a can to which the pacing lead is coupled. 

41. The GUI of claim 37, further comprising a field that indicates an 
impedance of the pacing lead. 

42. The GUI of claim 37, further comprising a field that indicates a 
temperature of a temperature sensor included with the pacing lead. 

43. The GUI of claim 37, further comprising a field that allows the user to 
specify a script file to run. 

44. The GUI of claim 37, further comprising an area for displaying in real 
time a signal detected by the pacing lead. 

45. The GUI of claim 37, wherein the area is configured to display an 
illustration of cardiac synchrony. 

46. The GUI of claim 45, 

wherein the illustration of cardiac synchrony is configured to trace a point in real 
time on an X-Y coordinate plane; 

wherein the X coordinate of the point indicates a value corresponding to the 
amplitude of a first response signal indicating cardiac contraction at a first location; and 
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wherein the Y coordinate of the point indicates a value corresponding to the 
amplitude of a second response signal indicating cardiac contraction at a second location. 

47. The GUI of claim 45, wherein the illustration of cardiac synchrony is 
configured to: 

5 trace a first point on a polar coordinate plane, wherein the radial coordinate of the 

first point indicates a value corresponding to the amplitude of a first signal indicating 
cardiac contraction at a first location, and wherein the angular coordinate of the first 
point is configured to change at a constant rate; 

trace a second point on the same polar coordinate plane, wherein the radial 
10 coordinate of the second point indicates a value corresponding to the amplitude of the 

second signal indicating cardiac contraction at a second location, and wherein the angular 
coordinate of the second point is configured to change at the same constant rate; and 
to highlight the area enclosed by both traces of the first point and the second 

point. 
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